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Abstract. Ill this paper, we consider some cubic near-Hamiltonian systems obtained from per- 
turbing the symmetric cubic Hamiltonian system with two symmetric singular points by cubic 
polynomials. First, following Han [2012] we develop a method to study the analytical property 
of the Melnikov function near the origin for near-Hamiltonian system having the origin as its 
elementary center or nilpotent center. Based on the method, a computationally efficient algo- 
rithm is established to systematically compute the coefficients of Melnikov function. Then, we 
consider the symmetric singular points and present the conditions for one of them to be elemen- 
tary center or nilpotent center. Under the condition for the singular point to be a center, we 
obtain the normal form of the Hamiltonian systems near the center. Moreover, perturbing the 
symmetric cubic Hamiltonian systems by cubic polynomials, we consider limit cycles bifurcating 
from the center using the algorithm to compute the coefficients of Melnikov function. Finally, 
perturbing the symmetric hamiltonian system by symmetric cubic polynomials, we consider the 
number of limit cycles near one of the symmetric centers of the symmetric near-Hamiltonian 
system, which is same to that of another center. 
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1 Introduction 

There have been many studies on bifurcations of limit cycles from elementary singular 
points, see Bautin [1952]-Yu & Han [2004]. In general, there are two types of such bi- 
furcations leading to limit cycles: either by perturbing a focus or by perturbing a center. 
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Especially, many authors study C°° systems of the form 



X = Hy + ep{x,y,S), y = -H^ + €q{x,y, 5), (1) 

where H{x, y), p{x, y, 6), q{x, y, 6) are C°° functions, e > is small and 6 E D C is a 
vector parameter with D compact. 
When e = 0, system ([1]) becomes 

X = Hy, y = -H^, (2) 

which is Hamiltonian system, and thus ([1]) is called a near-Hamiltonian system. We 
impose the following hypothesis 

(H) H-,{xo,yo) = Hy{xo,yo) = 

and denote A{xo,yo) = — — -^{xo,yo) and A = A(0,0). 

o[x, y) 

If the hypothesis (H) is satisfied, then system (j2]) has a critical point at the point 
{xo,yo). Suppose that the hypothesis (H) is satisfied at the origin. Then it is said that 
the origin is (i) a saddle if det{A) < 0; {ii) an elementary center if det{A) > 0; {Hi) a 
nilpotent critical point if det{A) = and A ^ 0. 

From now on, we always suppose that system ([2]) has a critical point at the origin. 
Without loss of generality, for det{A) > we may assume that the expansion of H at the 
origin is of the form 

H{x,y) = + + ^ h,,x'yJ^ ^ > 0, (3) 

i+j>3 

for det{A) = and A ^ we may assume 

H{x,y) = + ^^i^V' w > 0. (4) 

i+j>3 

We call ([3]) or (jlj) the normal form of Hamiltonian function of system (jj]) having the 
origin as its elementary center or nilpotent critical point respectively. Similarly, we call 
the corresponding Hamiltonian system or near-Hamiltonian system the normal form of 
system (E]) or ([T]) near the center. 

For the function H{x,y) of the normal form ([3]) or (jl]), the implicit function theorem 
implies that there exists a unique analytic function (p{x) = 0{x^) such that Hy{x, (p{x)) = 
for |x| small. Let 

H*{x) = H{x, ip{x)) = J2 hj^^^ hk^O, k>2. (5) 

Recently, Han, et al. [2010] gave a complete classification of nilpotent critical points for 
Hamiltonian systems (|2]) as follows. 

Theorem 1. The origin is (i) a cusp if k is odd; (ii) a saddle if k is even with hk < 0; 
and (Hi) a center if k is even with hk > 0. 



Definition 1. Let the Hamiltonian function H{x,y) be of the form or Then 
for system (j2]) the origin is called a cusp of order m if k = 2m + 1. It is called a nilpotent 
center of order m (a nilpotent saddle of order m, respectively) if k = 2m + 2 and /i^, > 
(if k = 2m + 2 and <0). 

Now suppose that the Hamiltonian system (j2]) has an elementary center or a nilpotent 
center of order p — 1 > 2) at the origin, namely the Hamiltonian function H{x, y) has 
the form (EI) or (jlj) and satisfies 

Hoix) = J2 ^i^"' H > 0> P > 1- (6) 

Then, the Hamiltonian system ([2]) has a family of periodic orbits, given by 

Lh:Hix,y) = h, /i G (0, /3) 

such that Lh approaches the origin as /i — )■ 0. 

Take h = Hq E (0, /3) and A^Hq) G Lho- Let / be a cross section of system ([21) passing 
through A{hQ). Then, for h near Hq the periodic orbit has a unique intersection point 
with I, denoted by A{h), i.e., A{h) = LkHl. Consider the positive orbit 'j{h,e,A) of 
system ([1]) starting from A{h). Let B{h,e,A) denote the first intersection point of the 
orbit with I. Then we have 

H{B)-H{A)= [ dH = €[M{h,5) + 0{e)]=£F{h,e,5), (7) 
Jab 

where 

M{h,6) = i {Hyq + H^p)dt 

' (8) 
= y {qdx - pdy) = jj {p^ + qy)dxdy. 

H<h 

The functions F{h, e, 5) and M{h, 6) in ([7]) are called bifurcation function and Melnikov 
function of system ([T]) near the origin respectively. The resulting map from A{h) to 
B{h,e,6) is called a Poincare map of system ([1]). Obviously, for small e system ([1]) has 
a limit cycle near the origin if and only if the function F{h, e, 6) has an isolated positive 
zero in h near h = 0. 

Based on the analytical property of the Melnikove function M{h, 6) and the number 
of limit cycles near the origin by the function, we have the following theorems (see Han, 
Jiang & Zhu [2008]). 

Theorem 2. Let system where H{x,y) has the form ([3]) or (jll), satisfy ([6]). Then, 
we have 

M{h,5) = h^^hj{5)hi. (9) 

i>o 

In particular, when p=l ( that is when system has an elementary center at the origin), 
we have 

M{h,6) = hJ2bj{S)h^- 

i>0 



Theorem 3. Under the condition of Theorem 2, if there exist k > 1, 6o G D such that 
bk{So) 7^ and 

O{0i, ■■ ■ ,dk) 

where 6 = {6i, ■ ■ ■ , 6m), m > k, then there exist a constant Eq > and a neighborhood V 
of the origin such that for some < |£:| < eo one? |5 — 5o| < (III) has k limit cycles in 
Vi. 

Theorem 4. Consider the near-Hamiltonian system where H{x,y,a) with a G R" 
satisfies ([6]) and the functions p and q are linear in S E R"^. Suppose there exist integer 
k > and Sq = (5io, ■ ■ ■ , 5mo) ^ R™; (^o ^ R" such that 

b^So, ao) = 0, J = 0, ■ ■ ■ , - 1, det " ' '\~'\ ao) ^ 0, (10) 

and 

h+j\{Si,-,5k)=ah+i,-,Sm,a) = Lj{Sk+i,--- ,6m)Aj{a), j = 0, ■■■ ,n, (11) 

where 

(Si, ■ ■■ ,Sk) = ^{Sk+i, ■ ■ ■ , 5 m, a) 
is the unique solution to bj = 0, j = 0, ■ ■ ■ ,k — 1 for a point a near aQ, 

Lj{Sk+i,o, ■ ■ ■ ,Smo) 7^ 0, j = 0, ■ ■ ■ ,n, 
Aj(ao) = 0, j = 0, - ■ ■ ,n - 1, A„,(ao) 7^ 0, 
and 

det ^^^;^---^\\ ao)^0. (13) 
d(ai, ■ ■ ■ ,a„) 

Then, for some {e,6,a) near (0,(5o,ao), system has k + n limit cycles near the origin. 

In light of the above theorems, a key step in studying the smaU-amphtude limit cycle 
bifurcations of system ([1]) is to find an efficient method to compute the coefficients bi. For 
p=l, the formulas for the first three coefficients bj{6), j = 0, 1,2 were obtained by Hou 
& Han [2006] by using the double integral in ([8]). Recently, Han, Yang & Yu [2009] and 
Han [20 12] have developed a new approach to prove Theorem 2 and have established a 
computationally efficient algorithm to systematically compute bj{6), j = 0, 1, 2, 3, ■ ■ ■ for 
p = 1 and p > 1 respectively. Following Han, Yang & Yu [2009] and Han[2012], we prove 
Theorem 2 and then also establish a computationally efficient algorithm to systematically 
compute bj{S), j = 0, 1, 2, 3, • ■ ■ for arbitrary positive integer p. 

In this paper, to illustrate the efficiency of the approach, we consider the following 
symmetric cubic Hamiltonian system 

(S): Hamiltonian system ([2]), with Hamiltonian function of the form 



(12) 



Hix,y)= ^ hijx'y\ 

2+1=2.4 



taking (0, ±1) as its critical points, 

We will consider the conditions for the symmetric points to be elementary centers or 
nilpotent centers. Under the condition for one of the symmetric points to be a center, 
perturbing system [S) by cubic polynomials one can obtain a cubic near-Hamiltonian 
system. Then, by suitable change of the variables (x, y) and t, we gave out the normal 
form of the above near-Hamiltonian system near the center (0, 1) (or (0, —1)). Obviously, 
the number of limit cycles of the near-Hamiltonian system near the center (0, 1) (or 
(0, —1)) is the same as that of the corresponding system of the normal form near the 
origin. As the result, we can study the number of limit cycles near the center (0, 1) 
(or (0,-1)) by studying the Melnikov of the normal form near the origin. Moreover, 
perturbing the symmetric hamiltonian system (S) by symmetric polynomials, one can 
obtain a symmetric near-Hamiltonian system. Then, if one found some limit cycles near 
one of the symmetric centers, then the system has the same number of limit cycles near 
the other center. Hence, it is easier to find some limit cycles of the system near both 
centers than to find the same number of limit cycles near one center. 

This paper is organized as follows. In the next section, we give a proof of Theorem 
2, and list an efficient algorithm for any positive integer p based on this proof as the 
appendix. In section 3, we study the critical points (0, ±1) of system (S) and list all 
conditions for one of the singular points to be an elementary center or a nilpotent center. 
Moreover, near the center we present the corresponding normal form of the Hamiltonian 
system (S). Then, in section 4, perturbing system {S) by cubic polynomials, we obtain 
the corresponding near-Hamiltonian system. Applying the theorems and the program to 
the normal form of the above near-Hamiltonian system, we estimate the number of limit 
cycles bifurcating from the center (0, 1) or (0, —1). In section 5, we discuss the number 
of limit cycles near both the symmetric centers of the symmetric cubic near-Hamiltonian 
systems obtained from perturbing system (5) by symmetric cubic polynomials. 

2 Proof of Theorem 2 

Recently, Han, Yang & Yu [2009] developed a new approach to prove Theorem 2 for p = 1 
and then established an algorithm based on its proof. The algorithm was implemented in 
the computer algebra system Maple. Earlier, Han, Jiang & Zhu [2008] proved Theorem 
2 for p > 1. Then following Han, Yang & Yu [2009], Han [2012] gave a new proof. In 
this paper, for convenience and as a preliminary, we repeat the proof of Theorem 2 for all 
positive integer p > 1 following Han [2012], and then establish an algorithm based on the 
proof and provide an implementation in computer algebra system Mathematica. 

To prove Theorem 2, we first introduce a change of variables to make the form of the 
Hamiltonian function simpler. By ([3]) and (jl]) using the implicit function theorem, one 
can show that there exists a unique C°° function (p{x) such that Hy{x,(p{x)) = for |x| 
small. Thus, we can write 




(14) 



By introducing a new variable v = y — f{x), system ([T]) can be rewritten as 

X = H*{x, v) + ep*{x, f , 5), 
V = -H*{x,v) + eq*{x,v,6), 

where 



(15) 



H*{x,v) = H{x,v + ip{x)), 

p*{x,v,6) = p{x,v + (p{x),6), (16) 
q*{x, V, 6) = q{x, V + (p{x), 6) — (p'{x)p*{x, v, 6). 
Noting that Hy{x, ip{x)) = 0, we have 

H*{x, v) = H*{x) + H*{xy+^ = H*{x) + v^H{x, v), (17) 
i>i 

where ^'(O, 0) = ^ > and 



H*{x) =H{x,^{x)) = Y,h 

1 d^+'H~' ^^^^ 

Let ([6]) be satisfied. Then it follows from (fT7|) that there exist a family of periodic orbits 
surrounding the origin defined by the equation H{x,y) = h or H*{x,v) = h for h > 
small. Let w = \/h — Hq{x) and suppose that the region surrounding by the closed curve 
H*{x,v) = h can be expressed as the form {{x,v)\x2{h) < x < xi{h), V2{x,w) < v < 
Vi{x,w)}. Then, from Han, Yang & Yu [2009] or Han [2012], we have 

Lemma 1. (?) The equation H*{x,v) = h has exactly two C°° solutions Vi{x,w) > and 
V2{x,w) < in V satisfying 

vi{x, w) = a/2w(1 + 0{\x, w\)), V2{x, w) = vi{x, —w). 

[a) Xi{h) > and X2{h) < are the two solutions to equation Hq{x) = h. 
Then, it is obvious that 

M{h, 6) = J) q*dx - p*dv 

JH*{x,v)=h 

{pI + ql)dxdv (19) 

H*<h 

xi{h) /•vi(x,w) 

/ iP*x + q*v)dvdx. 



Let 



q{x,v,S)= {pl + ql)dv = q*(x,v,S)-q*{x,0,S)+ pl{x,u,S)du. (20) 

In In 



Then, we have 



where 



q{x,v,5) = V hijX^v^ = '^^qj{x)v^ , (21) 
i+j>0 j>x 

1 
1 

(p, + g,)(x,¥p(x),0,(5) (22) 



i>0 



Further, we write 



By Lemma 1 we have 



vi{x,w) = ^^aj(x)w-^. (23) 



h-H*{x) = H*{x,v)- H*{x) = {vi{x,w)fJ2^j(^)lMx^'^)y~^ 

= aj{x)H*{x)w^ + {2ai{x)a2{x)H*{x) + al{x)H;{x))w''^+ 

[{al{x) + 2ai{x)a3{x))H*{x) + 3al{x)a2{x)H;{x) + af{x)H;{x)]w'^ + ■■■ . 
Equahng the coefficients of in the above identity, we have 

Then, it follows from fll9p that 

M[h,5) = I [q{x,Vi{x,w)) — q{x,V2{x,w))]dx. 

Jx2{h) 

By Lemma 1, the function q{x^vi) — q{x,V2) is odd in w. Thus, we can write 

q{x, vi) - q{x, V2) = ^ gj(x)w^^+\ (25) 
i>o 

and hence, 

/•xi(h) 



M{h,6) = J2 qjix)w^^^^dx. 

j>Q Jx2{h) 



'X2{h) 

In order to compute the above integral, we change the limits of integration. Let tplx) = 
j_ J- 
sgn(x)[ifo(x)]2p. Then, by (El) the function ip is C°° for small |x| with i>'{0) = h"^ > 0. 



Therefore, we may introduce the new variable u = ijj{x) to obtain 

1 

^ qj{u)w^^+^du 
= [qji^] + qji-u)W~^^du, 



1>0 



where w = \/h — u^p and 



It is easy to see that we can assume 



(26) 



q,{u) + q,{-u) = J2rijU^'- (27) 

j>0 

Then, 

M{h,6)= J2 ^vW, (28) 

i+j>0 

where 

1 1 

Iij{h)= / uV^'du= / u^'{h-u^P)Wh-u^Pdu. 
Jo Jo 

Lemma 2. Let 

1 

2i 



P^.= / v~{l-vy^l-v'^dv. (29) 
Jo 

Then 

Ii,[h) = f3ijh^^—, h>0, < Ai < 1- 
Proof. Introducing v = u^/h^, we have 

Iij{h) = h ^ p / VP {1- vyVl - v'^dv. 
Jo 

This ends the proof. 

Now by fl28l) and Lemma 2, we have 

M{h, 5) = h^ rijPijh-^ = h'^J2 (30) 

i+j>0 l>0 

where 

bi{S) = ''^^P^y (31) 

i+pj=l 

Finally based on fl2B|) . fl2B|) and fl27]) . we see that if system ([T]) is analytic, then the 
series ^ rijU^'^w'^^^^ is convergent for (m, w) near the origin. Then it follows that the 

series ^ |rjj|/i*"'"-' is convergent for some constant /i > 0, and hence by (15(1 M{h,6) is 

analytic in /i. 

This completes the proof of Theorem 2. 



3 The symmetric critical points of symmetric Hamil- 
tonian system [S) 

In this section, we study the symmetric critical points (0, ±1) of symmetric cubic Hamil- 
tonian system (S). Obviously, we have Hy{0, ±1) = Hx{0, ±1) = 0, i.e., 

±{2ho2 + 4/io4) = ±{hu + /iia) = 0. 

Thus, we can take 

ho4 = -ho2/2, hi3 = -/ill. (32) 

For the symmetry, we need only to consider the critical point (0, 1). By introducing a 
new variable v = y — 1 and then taking v as y, we have 

x = Hy{x,y), y = -H^{x, y), (33) 

where H{x,y) = H{x,y + 1). Truncating the constant term of the function H, we have 

H{x, y) = (/i2o + h22)x'^ - 2huxy - 2/io2?/^ + hix^ + 2h22X^y - Shuxy"^ 

1 (34) 
-2ho2y^ + h^ox^ + hsix^y + h22x'^y'^ - huxy-^ - -ho2y^. 

As a result, the point (0, 1) of system (S) corresponds to the origin of system (155]) . If 
(|32|) is satisfied, then system (|33|) has a critical point at the origin. Moreover, for system 
(l33l) we have 

/ -2/iii -4ho2 \ 
A=[ (35) 

\^ -2(/l20 + h22) 2hu J 

and 

det{A) = -Ahl^ - 8/io2(/i2o + h22). (36) 

Then, the origin is {i) a saddle if det{A) < 0; {ii) an elementary center if det{A) > 0; 
{in) a nilpotent critical point if det{A) = and A ^ 0. 

If A = det{A) > 0, which means that hQ2{h2o + /i22) 7^ 0, then system (133|1 has an 
elementary center at the origin. Let u = a/A > 0. Without loss of generality, we can let 
^20 + h22 = 1/2 and a; = 1 (otherwise, we can introduce a suitable rescaling of (x, y) and 
time t). Therefore, we obtain 

1 2 1 

h2Q = —h22 + - 5 hQ2 = —hii — - . 

Then, introducing the change of variables 

u = -y, V = x- 2huy, (37) 
and then taking {u,v) as {x,y), system (1551) becomes 



ho = -7;- 4(2/i22 + 2/111/131 - l)hl^, ho4 = hio, his = -hi - 8hih. 



40, 



where 

Hi{x,y) = -{x' + y^)+ ^^i^V, (39) 

3<j+i<4 

and 

^03 = ^31, = -2/i22 - Ghihi, hi = -^hi + 9>hih2 + I2hl^hi, 
1 
2 

h2 = h2 + Qhihi + 2Ahl^ho, hi = hi - 'ihih2 - l^hl^hi - ^2hliho, 
13 

^40 = ^ - TT^ii + 4/i?i/i22 + Shl^hi + IQhj^ho- 
o Z 

Obviously, if the symmetric Hamiltonian system (S) has an elementary center at (0, 1), 
then the normal form of system [S) near (0, 1) should be system ( l38l) . For the symmetry, 
system {S) has an elementary center at (0, —1) and is of the same normal form fl38|) at 
the same time. 

Now, suppose that system f l33|) has a nilpotent critical point at the origin. In other 
words, we have 

- Ah,, - 8h2iho + h2) = 0, + h^2 + iho + h2y ^ 0. (40) 
If h2 = 0, then hi = and ho + ^22 7^ 0. Thus, system (!33|) becomes 
X = x^(2/i22 - hix - 2/1222/) 

y = ~x[2{ho + h2) + ^hix + 4/1222/ + ^hox"^ + "ihixy + 2h2y% 

which means that all the points on the line x = are singular. Then the origin is not an 
isolated singular point. 

Therefore, we can assume that h2 7^ 0. Without loss of generality, suppose h2 = - 
(otherwise, we can introduce suitable rescaling of {x,y)). Then, from (HOl) we can take 

h2 = -/^n - ho- (41) 

Introducing new variables as 

1 , 

u = --x, v = hix + y, (42) 
and then taking {u,v) as {x,y), system (l33l) becomes 

where ^ 

H2{x,y) = -y^+ hj^y^ (44) 

3<i+i<4 



^30 = 8hnh2o + 4/2,31, h2i = -2/1^ + 4/i2o, hi2 = = 0, 

ho3 = ^, = 2/iii(/in + 4/120) + 8(/iii/i3i - h^o), (45) 

/i31 = 8/lii/l20 + 4/I31, /;,22 = + 2/l20, /i04 = ^ • 

o 

Using the mathematica code we compute 

62 = -4/i2o + 2hj^ , 

63 = -4(2/;,ii/;,2o + /^si) , 

64 = -2(/i2, - 2M' , 

(46) 

65 = 16(/iii - 2/i2o)(2/iii/i2o + ^31) , 

66 = 4(/ifi - 6/iti/i2o - l2hl,hlo - 8/1^0 - 24/in/i2o/i3i - Qh^) , 



and 



^3 

/l6 



4(2/in/i2o + hsi) , 

8hn{2hnh2o + h^i) - Khlo + ^4o) 

8(/l2,- 2/102) (2/lll/l20 + /i3l) , 

-8(2/in/i2o + /i3i)' , 



(47) 



Introduce the following conditions: 

(A) /i3i ^ -2huh2o; 

(B) hsi = -2/iii/i2o, ^io + ^40 > 0; 

(C) /i3i = -2/iii/i2o, /iio + ^40 < 0; 

(D) h-ii = -2/iii/i2o, = -h^Q. 

Then, we can prove the following statement. 

Theorem 5. For system . the origin is an isolated nilpotent singular point if and 
only if one of conditions (A), {B), (C) listed above holds. Further, the origin is 

(i) a cusp of order 1 if and only if {A) holds, 

(ii) a nilpotent saddle of order 1 if and only if (B) holds, 
(Hi) a nilpotent center of order 1 if and only if (C) holds. 

Proof. By Definition 1, the origin is a cusp of order 1 under {A). If [B) or (C) holds, 
then /is = 0, < or = 0, /;.4 > 0, and the origin is a nilpotent saddle or center of 
order 1. If {D) holds, then we have 

^2(x, y) = l{y^ + 2y~ Ahj.x^ + 8/l2ox')^ 



which means that all points on the curve 

+ 2y- Ahl^x^ + 8h2oX^ = 

are singular. This ends the proof. 

Assume that the origin is a nilpotent center of order 1 of system ( 1431) . i.e., h^i = 
— 2/iii/i2o and h^Q + h^o < 0. Then system ( H3|) becomes 

where 11 1 

y) = 2 2/' + 2^^'2/ + 2 2/' + Bx^ + ^^'2/' + g?/' (49) 

with 

A = 2h2o - B = 2h\^ - Shl^h^o - 8/140 

satisfying B > 2 . 

Similarly, if the symmetric Hamiltonian system (S*) has a nilpotent center of order 1 
at (0, 1), then the normal form of system (S") near (0, 1) should be system ( HHj) . For the 
symmetry, system (S*) has a nilpotent center of order 1 at (0,-1) and is of the same 
normal form fHHll. 



4 Bifurcation from the centers of systems [S] 

From section 3, we know that the symmetric Hamiltonian system (S") can take (0, ±1) as 
elementary centers or nilpotent centers of order 1. Under the conditions for one of the 
symmetric critical points to be a center of system {S\ we perturb it by cubic polynomials 
to obtain a near-Hamiltonian system and then study limit cycles bifurcating from the 
center (0,1) or (0,-1). 

Before starting to consider limit cycles bifurcating from a center of the near-Hamiltonian 
systems, we present some results on how to calculating the coefficients of the Melnikov 
function near the center. Suppose that system (|2]) has a center at the point (xq, 2/0) and 
H{xo,yo) = ho- Then, there exist a family of periodic orbits given by 

Lh : H{x, y) = h, he {ho, /iq + /?), /? > 0. 

To transform the linear part of the Hamiltonian system at a center into a normal form, 
one often introduces suitable linear change of variables or rescaling of the time. This 
procedure may cause a change of the Melnikov function. More precisely, making a linear 
change of the form 

u = a{x- xo) + b{y - yo), v = c{x - xq) + d{y - yo), 

and rescaling time t = k t, where D = ad — be ^ 0, system ([1]) becomes 



where ^ 

H{u,v) = —[H{x,y) - ho], 

1 

p{u, v,S) = - [a p{x, y,S) + b q{x, y, 6)], 
q{u, v,6) = ^ [c p{x, y,6) + d q{x, y, 6)], 

X = xo + -^{du - bv), y = yo + ^("Cm + av). 
Let M denote the Melnikov function of system fl50|) at the origin, which is given by 

M{h,6)=(f) qdu-pdv\e=o, he{0,f3'). 

JH{u,v)=h 

Then, it is easy to see that 

x,y)=ho+kh/Dl^^^ cp)dx + {bq dp)dy]\e=o 
= sgnik)j §^^^^y^^^^^^^^^^[qdx - pdy] \ ,=o 

= ^^Miho + ^h,6), /iG(0,/3'), 

where M{h,6) is the Melnikov function of system ([1]) near the center {xo,yo) of the 
corresponding Hamiltonian system ([2]). Therefore, for < h < ho + f3 we have 

M{h,5) = ^M{j{h-ho,5)). (51) 



By Theorem 2, we have 

Theorem 6. Let system (I3U|) . where H has the form ^ or (jl]), satisfy Then, for 
ho < h < ho + (3 we have 

p+i 



M{h, 6) = {h- ho)—J2bA^)ih - ho)K (52) 



where 



m = sgn{k)[-) b,i6). 



Especially, for p=l, namely system ([2]) has an elementary center at the origin, then, 
for ho < h < ho + (3 we have 

M{h,S) = {h-ho)Y,b,i^)ih-hoy, 

j>0 

where 

m = sgn{k)[^yb,{5). 



As we known, D ^ Q and /c 7^ 0. Then, by Theorem 3 and Theorem 4, to determine 
the number of hmit cycles bifurcating from the center (xq, yo) of system ([T]), we need only 
to consider the Melnikov function M of near-Hamiltonian system ( l5Up . 

Theorem 9 told us a method to study limit cycles bifurcating from a center (xq, Uq) of 
the near-Hamiltonian system ([1]). First, by introducing suitable linear change of variables 
or rescaling the time, we can transform system ([1]) into a normal form (!50|) near the center 
(xo, yo). Then, using Theorem 2 after executing the Mathematica program, we can obtain 
the first several coefficients of the Melnikov function of system f l50|) at the origin. Finally, 
by Theorems 3 and 4 we can find the number of limit cycles of near-Hamiltonian system 
(l50l) at the origin, which is the same to that of the near-Hamiltonian system ([T]) at the 
center (xo,i/o)- 

Suppose that for systems (II]) and (l50l) . 



2 2 

Px + qy= ^ CijX^y^ Pu + Qv = ^ CijUV. 
i+j=0 i+j=0 

Then, using the Mathematica program we can compute 

Coo = ^[coo + xq cio + yo coi + C20 + xoyo cn + y^ C02], 
1 

Cio = ('^01 + xo Cn + 2yo C02) + d (cio + Vo Cn + 2xo C20)], 

coi = ^['^ ('^01 + ^0 Cn + 2yo C02) - b (cio + yo cn + 2xo C20)], 

C20 = ^jj2 t*^^ '^02 — C(i Cn -|- (i^ C20], 

Cn = ^^t^^*^*^ ^02 + {be + ad) Cn - 2bd C20], 

C02 = ^^[a^ C02 - a6 Cn + C2o]. 

Now, we turn to considering the near-Hamiltonian system obtained by perturbing sys- 
tem (S) with cubic polynomials as follows: 



(53) 



where 



X = Hy + ep{x,y), y = -H^ + eq{x,y), (54) 
Hix, y)= ^ hijx'y^ 



i+j=2,4 



p{x,y) = "y^V, ?(a;,i/) = "Y b, 

i+j=l i+j=l 



tjX y'^ 



and 



Px + qy = Coo + CloX + Coiy + C2qX^ + CxxXy + C2o2/^. 



From section 3, if A = — 4/iJ^ — 8/102(^20 + ^22) > 0, which means that /;,02 7^ and 
/i2o + ^22 7^ 0, then there is an elementary center of Hamiltonian system iS) at the point 



(0, 1). For convenience, we let /i2o + ^22 = 1/2 and A = 1. Then, introducing a change of 
variables (x, y) as 

u = -{y - 1), V = x- 2hn{y - 1), 

and then taking {u,v) as {x,y), one can obtain the normal form of system (15^ near the 
center (0, 1) as 

du dHi . , ^ , ^ dv dHi . , ~ , , /^^n 
'dt^~dv~ ^'^ = — -q;^{u,v) + £qi{u,v), (55) 

where i^i is the same as that in section 3 and 

2 



dpi , <9gi ^ 



where 



Cloo — Coo + Coi + C02, Clio — — 2/liiCio — Cqi — 2Co2 — 2/liiCii, cloi — Cio + Cii, 
CI20 = C02 + 2/liiCii + Clii = -Cii - 4/liiC20, clo2 = C20. 

Obviously, cIqo, clio, cIqi, CI20, cln, clo2 are all linear functions of Cqo, Cio, Cqi, C20, Cn, C02 
and are independent. Therefore, we can denote cljj as for 0<i+j<2,i>0, j>0. 

For system (!55l) . there are too many parameters in the corresponding Hamiltonian 
function. Hence, the coefficients of the Melnikov function will be very complicated and 
it is very hard to give out the accurate number of limit cycles of system fl55l) near the 
origin. However, suppose that there are altogether m parameters in (5, a), by Theorem 3 
or Theorem 4 we know that one can find at most m — 1 limit cycles near the origin. As 
the result, one can find at most 9 limit cycles of system (l55ll near the origin. Here, we 
consider an simple case of system (l55l) to illustrate that system fl55|) can have at least 6 
limit cycles near the origin. As an especial case, we assume that hu = /i22 = 0, /131 = 1 
in system ( l55i) . Namely, we consider the following system, 

ii = Hy + ep{u, v), i) = —H^ + eq{u, v), (56) 

where p, q are the same as before and 

H{u, v) = -{u'^ + v"^) - -M^ + - - uv^ + rv^ 
2 2 8 

with r = h^Q. 

Theorem 7. For system (l56l) . introduce 

^ = (coO) Cio, Cqi, C02, Cii, C20), c = (r), 

and 

"0 — [Cqo, CiQ, Cqi, Cq2, Cii, C2q) , CTq — {r ), 



where 5q and ctq satisfy 

3 

/(r*) = 0, Cqo = 0, = ~2^*o ~'~ "^^01 ~ ^20' '^oi 7^ 0' 
= (12r* - ^)cto + (40r* + 3)c*i + (8r* - ^)c*o, 
= 17680r* - 43347 [(^^^^(r*)^ - 17328r* - mQ)cl, + (+11107 - 5680r*)c^o] , 



6544683 + 11682846r* + 30378720(r*)2 - 40278528(r*)3 - 21471232(r 



*\4 



''^^ 2(2770317 - 11417834r* + 11442688(r*)2 + 442880(r*)3) 

and 

f{r) = 780900831-459924741r-2093583104r2-1597844992r^-108363776r^+1381957632r^ 

Then, for some [e, 6, a) near (0, 6q, ctq) system fl56|) has 6 limit cycles near the origin. 

Proof. Let M{h,6) be the Melnikov function of system {nH)i near the origin. Then, by 
Theorem 2, we have 

M{h,5) = hJ2bj{S)h^. 

i>o 

Executing the Mathematica program yields 

bo = 2 7rcoo. 

Letting cqo = 0, under which Bq = 0, we obtain the following expressions for 6j's. 

3 

&i|6o=o = 7r[c2o + Co2 3coi + 2*^10]' 



^2 1 60=0 



^ [(25 - 12r)cio + 140(-3 + 2r)coi + (25 - 4r)c2o - 4cn + 8(18 - 5r)co2] ^ 



bi\bo=0 — hi^ilClO + ^i2Coi + 5j3C20 + '^i4Cll + 5j5Co2); i — 3, 4, 5, 6 

(57) 



where 



5n , _ 2l7r 2l7r _ 337r 

, La — , I's — 1 tfi — ■ 

64' 512' 2048 8193' 



and 

= 939 _ 210r - 84r2, = 462(39 - 52r + 12r^), 833 = -2(333 - 270r + 28r'^), 
534 = -2(123 + 28r), ^35 = -8(749 - 696r + 63r^), 

^41 = 85233 - 85896r + 6584r2 + lOSGr^, 5^2 = 3432(-331 + 646r - 340r2 + AOr^), 
Si3 = 21278 - 34416r + 13968r2 - 704r3, S^^ = -8(-3009 + 1706r + 132r2), 

= 64(-5828 + 9101r - 3222r2 + 143r3), 
551 = 13788577 - 23251746r + 9251880r2 - 382800r3 - 34320r^ 

= 92378(1821 - 4664r + 3864r2 - 1120r3 + 80r^), 
^53 = -2(1156991 - 2579166r + 1845528r2 - 391600r3 + 11440r^), 
554 = -2(2028023 - 2589736r + 582120r2 + 22880r3), 

= -32(1717743 - 3736392r + 2422830r2 - 462176r3 + 12155r^), 
561 = -1454668039 + 3393291252r - 2432266980r2 + 541324640r3 - 13999440r^ - 806208r^ 
5q2 = 386308(-44055 + 139650r - 157944r2 + 75600r3 - 14000r^ + 672r5), 
5g3 = 2(97965305 - 271364748r + 269902620r2 - 108076960r3 + 14155440r4 - 268736r5), 
^64 = 4(109528617 - 212220190r + 113651592r2 - 14093040r3 - 335920r^), 
^65 = 128(43337762 - 120692295r + 114411600r2 - 42231670r3 + 5071950r^ - 88179r^). 

Let 

6, = 6,|6o=o, J = 1,2,3,4,5,6. (58) 
Noticing ([57]), we can solve C02, Cn from 61 = 62 = as follows: 

3 191 47 

C02 = -2^10 + 3coi - C20, Cu = (12r ^)cio + (40r + 3)coi + (8r - y)c2o- (59) 

Substituting fl59|) into 63, 6465 and be results in 

b3= — r(17680r - 43347)cio + 16(87 + 2166r - 224r2)coi + 2(-11107 + 5680r)c2o] , 
128 

h = ^(^iiCio + S'i^coi + 5-3C2o), i = 4, 5, 6 

(60) 



where k same as before and 

5^1 = -2(1794159 - 1900148r + 429248r2), 6'^^ = 64(1725 + 43488r - 34292r2 + 2112r3), 

= -4(458703 - 530676r + 138304r2), 
5'^^ = 579832875 - 997181804r + 486587688r2 - 68706176r^ 

= 16(-1106391 - 27357598r + 41005384r2 - 13225168r'^ + 549120r^), 

= -2(-147970875 + 268572172r - 142634920r2 + 22244992r3), 
S'^^ = -4(15347742095 - 36178820899r + 28146658686r2 - 8398426168r=^ + 800059520r^), 
(5^2 = 32(58577817+ 1411250460r - 3104524228r2 + 1937691360r=^ - 357475040r^ + 10749440r5), 
5^3 = -8(3911748231 - 9589874467r + 7894793822r2 - 2537774136r3 + 260074880r^). 

For r near r* we have 17680r - 43347 0. Thus, 63 = if and only if 

cio = i7680r- 43347 t^^^^^""' ~ ~ ^^^^"""^ + (+11107 - 5680r)c2o] , 

which results in 

64 = — r [(6544683 + 11682846r + 30378720r2 - 40278528r3 

32(-43347+ 17680r) 

-21471232r^)coi - 2(2770317 - 11417834r + 11442688r2 + 442880r3)c2o] ■ 

For r near r* we have 2770317 - 11417834r + 11442688r2 + 442880r3 7^ 0. Thus, 64 = 
if and only if 

_ 6544683 + 11682846r + 30378720r2 - 40278528r3 - 21471232r^ 
~ 2(2770317 - 11417834r + 11442688r2 + 442880r3) 

which results in 

= AjL, i = 0, 1, 

where 

37r(4r^ + 4r - 7)coi 



32(2770317 - 11417834r + 11442688r2 + 442880r3) ' 

Ao = -14/(r), 

Ai = ll(-77526899979 + 74770117515r + 170378796326r2 + 50680931328r3 
-71020059648r^ - 138129047552r^ + 61673963520r6). 

Obviously, there is at least one real solution to equation /(r) = and denote it by r*. 
Then, for r near r* and Cqi near Cq^ we have L ^ 0. Thus, 65 = if and only if Aq = 0. 
When r = r*, we have 65 = but Ai 7^ 0, which results in ^ 0. Hence, it is clear that 
the conclusion follows from Theorem 4. 



The proof is complete. 



Therefore, if system (S) has an elementary center at (0, 1) (or (0, —1)), then one can 
find 6 to 9 limit cycles of system (15^ near the center (0, 1) (or (0, —1)). 

If A = —4:hl^ — 8/102(^20 + ^22) = 0, /i3i = — 2/iii/i2o and /i^o + ^40 < 0, then there is a 
nilpotent center of order 1 of Hamiltonian system (S) at the point (0, 1). For convenience, 
let ho2 = 1/2. Then, we can take 

h22 = —h\i — /l20- 

Introducing a change of variables (x, y) as 

u = ~{y - I), V = hiix + y-l, 

one can obtain the normal form of system (IMll near the center (0, 1) as 
du dH2 _ , dv dH2 

- = ^{u,v) + ep2{u,v),- = -^iu,v)+eq2iu,v), (61) 
where H2 is the same as that in section 3, namely, 

H2{u,v) = lv^ + 2{2h2o-hl,)u^v + lv^ + {2hl^-8hl^h2o-8h4oW + {2h2o-hl^yv^ + lv\ 
2 2 o 

and 



iU V- 



^ + ^ = Y c2, 

i+j=0 

where 

c2oo = Coo + coi + C02, c2io = -2cio + 2huCoi - 2cu + 4/iiiCo2, c2oi = cqi + 2co2, 
c22o = 4c2o - 4:hiicu + 4:hl^co2, c2u = -2cn + 4/iiiCo2, c2o2 = C02. 

Obviously, c2oo, c2io, c2oi, c22o, c2ii, c2o2 are all linear functions of (coo, cio, cqi, C20, cn, C02) 
and are independent. Therefore, we can denote c2jj as Cij for < i + j < 2, i > 0, j > 0. 



Theorem 8. For system (!6T|) . we introduce 

6 = (coo, cio, coi, C02, cii, C20), o- = {A, B), 

and 

Oq — \Cqq,Ciq,Cq-^^,Cq2,C^i,C2q), <7q — [/i , n ), 

where Sq and ao satisfy 

A* = 0, B*> 0, c*o = 0, c*2o = 0, CS2 = ^c*i, c*i ^ 0. 
Then, for some {e, 6, a) near (0, 60, Cq) the system (|6T|) has 4 limit cycles near the origin. 



Proof. Let M{h, S) be the Melnikov function of system (!6T!) near the origin. Then, by 
Theorem 2, we have 

M{h,6) = h"^^bj{5)yl^. 

j>o 

Executing the Mathematica program yields 



2r(5/4)V2^ 
° ~ r(7/4)(5- 2^2)1/4^00- 

Letting cqo = 0, under which = 0, we obtain the following expressions for biS. 



r(9/4)(5- 2^2)3/4 

r(i/4)y^ 

2r(ll/4)(fi- 2^2)5/4 



^3ko=o = r(13/i[ff^2)r/4 [2^(^' - ^)(5coi - 4co2) + (^2 + B)c,o], 

= ~ i6r(iV4)(^^^A2)^/4 [(^-^' - 6^'^ + 3^')^^^°^ - s^'^^) + 2^^^' - 3^)^^°]- 

(62) 

Let 

bj = bj\bo=o, j = 1,2,3. (63) 

Noticing ( 1^ . we can obtain C20 = 2Acoi by solving &i = 0. Substituting C20 = 2Acoi into 
62, 63 and 64 results in 



= r(ll/4)Ti^-S^)V4 [(5-4^ - 35)coi - 2(A^ - 5)co2], 
63 = ^4^jMh/2^ 2 _ 2B)coi - 2(^2 - 5)co2], 

r(13/4)(5 - 2^2)7/4 ^ ^ 02J, 

h = 15r(l/4)v^ _ ^ _ ^.^4 _ ^25 + 3E2)C02)]. 

16r(15/4)(fi- 2^2)9/4 ^ ^ ' "^'^ 

For (A, 5) near {A*, B*), we have ^2 - 5 < 0. Then 62 = if and only if 



which yields 



_ 5^2-35 

- 2(^2 _ 5)^01, 



~ _ 3Ar(3/4)v^ 
^ ~ r(13/4)(5- 2^2)3/4^01' 

b- 6r(9/4)v/2^ 942R,R2v 

- "r(15/4)(A2-i?)(5- 2^2)5/4(3^ - 2A 5 + 5 )coi. 



Thus, it follows that when 63 = we have A = or cqi = 0. However, cqi is a common 
factor of 63 and 64. Taking A = A* = 0, we have 

~ _ 6r(9/4)v^ 

r(15/4)i?V4^oi- 

Hence, it is clear that the conclusion follows from Theorem 4. 
The proof is complete. 

Therefore, if system (S) has a nilpotent center of order 1 at (0, 1) (or (0, —1)), then one 
can find 4 limit cycles of system fl54p near the center (0, 1) (or (0, —1)). 



5 Bifurcation in the symmetric near-Hamiltonian sys- 
tems 

In this section, we will consider limit cycles bifurcating from the symmetric centers of 
some symmetric near-Hamiltonian systems. We said a near-hamiltonian system ([1]) is 
symmetric, which means that not only the corresponding Hamiltonian function is sym- 
metric, but also the perturbing terms p{x,y) and q{x,y) are also symmetric with respect 
to the origin. 

Now, we consider the following symmetric near-Hamiltonian system: 

X = Hy + ep{x,y), y = -H^ + eq{x,y), (64) 

where 

H{,x,y)= hijx'y\ 

i+j=2,4 

p{.x,y)= ^ aijx'y\ q{x,y) = ^ hjx'y^ 

i+j=l,3 i+i=l,3 

and 

Px + qy = Coo + C2ox'^ + CuXy + CaoZ/^- 

Obviously, system flM|) is obtained by perturbing symmetric Hamiltonian system {S) 
with symmetric cubic polynomials and is symmetric with respect to the origin. If there 
are m limit cycles bifurcating from the center (xo, ?/o) (iiot the origin) of the corresponding 
Hamiltonian system, then it is same to the symmetric center (— xo, —yo)- Namely, there 
are 2m limit cycles of the symmetric near-Hamiltonian system. 

From section 3, if A = — 4/;,^^ — 8/102(^20 + ^22) > 0, which means that 7^ and 
^20 + ^22 7^ 0, then there is an elementary center of Hamiltonian system {S) at the point 
(0, 1). For convenience, we let + = 1/2 and A = 1. Then, introducing a change of 
variables (x, y) as 

u = -iy-l), v = X- 2hn{y - 1), 



one can obtain the normal form of system near the center (0, 1) as 

du dHi . , , , ^ dv dHi , . ^ , , 

'dt^~dv~ ^'^ = — -q;^{u,v) + eqi{u,v), (65) 

where Hi is the same as that in section 3 and 

dpi , dqi A . . 
ou ov ^-^ 

i+j=0 

where 

Cloo = Coo + C02, Clio = -2co2 - 2/liiCii, Cloi = Cii, 

Cl20 = Co2 + 2/liiCii + 4:hliC20, Clii = -Cii - 4/liiC20, clo2 = C20 

For simphfication, we study system f l65|) in two cases /in = and /in = 1 separately. 
First, let /ill = 0. Then for system (l65l) we have 

11 1 
ifi(M, t;) = -{u^ + v^) - -M^ - 2h22uv^ + hiv"^ + -u'^ + h22U^v^ - hsiuv^ + h^ov'^ 
2 2 o 

and 

cloo = Coo + Co2; Clio = — 2Co2, cloi = Cii, 
CI20 = C02, Clii = — Cii, clo2 = C20- 

Therefore, we have the following theorem 
Theorem 9. Suppose system (1651) satisfies hn = and 

ho2{h20 + h22) < 0. 

Introduce 

^ = (coO) C02, Cii, C20), cr = (/isi, h22, h^o), 

and 

^0 = (Coo; Cq2, Cii, C20), o"o = {h^i, h22, h^o), 

where Sq and ctq satisfy 

Cqq = — Cq2, C20 = 3/^3^0^]^ + (2 + 4/l22)Co2, ^ 0, 

^3l[20(/^22)'- 8/^22 -10/^:0-l] „. 

°2 -2 + 20(/i^i)2 + 8/1^2 - 4/ilo 



Then, for some {e, 6, a) near (0, Sq, o"o) the system f l65|) /ias 5 /zmzt cycles near the origin. 

Proof. Let M(/i, 5) be the Melnikov function of system (jHSD near the origin. Then, by 
Theorem 2, we have 

M{h,6) = hJ2bj{S)h^- 

i>0 



Executing the Mathematica program yields 

60 = 2 7r (coo + C02). 

Letting cqo = — C02, which imphes 60 = 0, we obtain the following expressions for 6j's: 
&i|bo=o = 7r[c2o - 3/i3iCii - (2 + 4/122)002], 

&2|bo=0 = ^ [(1 + 4/^22 + 20/^22 + 35/^2^ - 10/140)C20 

+ (/i3i + 20/122/^31 - 140/1^2/^31 - 105/iii + 70/i3i/i4o)cii (66) 
+ (-10 + lO/i^i - 4/122(3 + 6/i22 + 20/1^2 + 35/iii - 10/l40) + 4/i4o)co2] , 

6i|feo=0 = Li{6iiC20 + <5j2Cll + (5i3Co2)5 = 3, 4, 5, 



where 



and 



T ^ T ^ T 21 

L-x = vr, La = vr, Lr, = vr, 

32 ' 64 ' 1024 ' 



3, = -7 - 8/l22(3 + /l22(9 + 14/l22(2 + 9/l22))) " 3003/i^i 

+12/140 + 112/122(1 + 9/i22)/^4o - 252/^^^ 

+ 14/lii(l + 36(1 - Il/l22)/i22 + 198/i4o) 

32 = /i3i(-3 - 24/i22 + 28/i4o + 21(-96/i^2 + 528/1^2 + 429/i^i, 

+2/^2^(3 + 44/122 (-3 + 13/122) - 286/l40) + 48/^22/^40 

+132/ilo-8/ii2(l + 66/i4o))), 

-2(-63 - 2/122(35 + 4/122(15 + 2/122(15 + 7/122(5 + 18/122)))) 

+231(3 - 26/l22)/^3l + 12/l40 + 24/122(3 + 14/122(1 + 6/l22))/i40 

-28(1 + 18/i22)/i4o - 6/iii(-l + 14/i22(-l + 6/i22(-3 + 22/122) 
-66/i4o)+42/i4o)), 



■^33 



Hi 



Let 

6, = 6,|6o=o, J = 1,2,3,4,5. (67) 
Noticing (l66ll . we can solve C20 from 61 = as 

C20 = 3/i3iCn + (2 + 4/i22)co2 (68) 

Substituting (l68l) into 62 results in 

62 = 27r [2{10hl^ + 4/1^2 - 2/i4o - l)co2 + /i3i (10/^40 + 8/i22 - 20/1^2 + l)cn] • 

For a near cr* we have lOh^i + 4/122 ~ 2/140 — 1 7^ 0. Thus, we see that 62 = if and only if 

^ /l3l(20/l22- 8/122 -10/^40-1) , . 

""'^ -2 + 20/i2 + 8/100 - 4/i.n ^ ^ 



Substituting (!68|) and (!69|) into 63, 64 and 65 results in 

h+^ = ^iL, 2 = 0,1,2, (70) 

where 

^ /i3i(l - + lQhlh22 - 8hl^ + 16/1^2 + 4/^40 - 16/1^2/^40 + 4/1I0) 

16(10/iii + 4/^22 - 2/i4o - 1) ''''''' 



and 



Ao = 80(1 + 12/i22 - 28/1^2 + 21/^31 + 14/i4o), 

Ai = 168(3 + 32/i22 - 40/;,|2 + 192/i32 - 528/i^2 + 36/iii 

+264/122/^31 + 429/i|i + 28/i4o - 96/^22/^40 

+528/ii2/i4o - 132/i2g), 
A2 = 21(143 + 1452/122 - 1340/1^2 + 9632/112 - 16368/1^2 + 45760/i^2 



-137280/1^2 + 1509/^2^ + 10648/122/^31 - UlGOhl^h 



2 ^,2 
31 



+205920/1^2/1^^ - 194480/i|2/i3i + 11869/1^1 + 48620/^22/^31 
+184756/i|2/i3i + 138567/1^1 + 1158/i4o - 4528/^22/^40 
+18832/i22/;,4o - 45760/1^2^40 + 205920/i^2^40 
+2860/i^i/i4o - 102960/122/^31^40 + 194480/i^2^ii^40 
-92378/ii/i4o - 5324/^2 + 11440/122/?-^ - 



31/1.40 ^o^-*/t4o -r -L-L«w,t22"'40 

^2/^40 - 48620/i^i/i^o + 17160/i^o) 



For Cii near cli and a near a* we have L 7^ 0. Thus, 63 = if and only if Aq = 0, which 
yields 



/i4o = -^(1 + 12/^22 - 28/1^2 + 21/^2 J 



In this case, we have 
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Ai = y [4(1 - 2/122)' + 420(-l + 2/z22)/z3i + 1617/i^i] . 
Then, Ai = if and only if /i22 = fii^zi) or /i22 = f2{h^i), where 

/,(/^3i) = \- ^^hi, - {-iruv?>hi,, z = 1,2. 

Thus, it follows that when /122 = fiih^i) we have 

A2 = 688128(54+ (-l)*3lV3)/i^i, i = 1,2. 

Hence, it is clear that the conclusion follows from Theorem 4. 
The proof is complete. 



Now we turn to the case hu = 1. Similarly, we can prove the following theorem 
Theorem 9'. Suppose system ( l65l) satisfies hn = 1 and 

-4 - 8/io2(/i2o + ^22) > 0. 

Introduce 

6 = (coo, C02, Cii, C20), Cr = (/I31, /l22, /i4o), 

and 

Sq = (Cqo? Cq2, Cii, C20), O"0 = (^31' ^22) ^40)5 

where 6^ and ctq satisfy 



hlo = - f (^22)^ + |^/^;2 - 12/^22/^31 + f /^31 - ?(/^^3)^ 

1 



'^00 ~ '^02' '^20 



5 



[(52/1*2 + 60/1*1 - 22)c*2 + (75/1*1 + GOh*,, - 2Q)cl,] , 0^2 ^ 0, 



— 2L*Cq2 



-11 



;-8 + 20h*22 + 25/i^i)(-263 + 20(42 - 25/1^2)^^22 + 500/i^i + 1250/i^o) 
wzi/i 

L* = 979 + 6250(/i^i)2 - 50/i^i(101 + 4/i*2(-72 + 25/1^2) - 250/i^o) - 4050/;,^o 
+4/i;2(-1366 + 5 * (457 - 200h*2)h;2 + 2500/i*o). 

Then, for some {e, 6, a) near (0, 60, ctq) the system (I65p /ias 5 Hmz^ cycles near the origin. 

Therefore, system fl64|) has 5 limit cycles near the elementary center (0,1). For the 
symmetry, system fl64p also has 5 limit cycles near the center (0,-1). Namely, there 
are altogether 10 limit cycles of system fl64|) bifurcating from both symmetric elementary 
centers (0, ±1). 

If A = —4:hl^ — 9>hQ2{h2Q + /i22) = 0, h^i = — 2/111/120 and H^q + h^Q < 0, then there is a 
nilpotent center of order 1 of Hamiltonian system [S) at the point (0, 1). For convenience, 
let /io2 = 1/2- Then, we can take 

h22 = ~h\i — /i20- 

Introducing a change of variables (x, y) as 

u = ~{y -1), V = hiix + y-l, 

one can obtain the normal form of system flM|) near the center (0, 1) as 

du dH2 , . ^ , . dv dH2 / n . / n /_-, n 

-TT = ^) + ^P2{u, v), -7- = — -K-KU, v) + eq2{u, V), (71) 



where H2 is the same as that in section 3, namely, 

H2{u,v) = ^l;2 + (4/,2o-2/^?l)n2^; + l^;3 + 2(/^4^-4/l?l/l2o-4/^4o)w' + (2/^2o-/^n)«V + ^^;^ 

Z Z o 

and 

2 

dp2 , dq2 



aw aw ^-^ 

where 

c2oo = Coo + C02, c2io = — 2cii + 4/iiiCo2, c2oi = 2co2, 

c22o = 4c2o - 4/iiiCii + 4/ifiCo2, c2ii = -2cii + 4/iiiCo2, c2o2 = C02. 

Theorem 10. Suppose system ( I7T]) satisfies 

-4:hl^ - 8/l02(/i20 + ^22) = 0, /l3i = -2/lii/l20, /i20 + ^40 < 0. 
(5 = (coo, C02, Cii, C20), CI = (/ill, /i20, /i4o), 

and 

<^0 = (CoO) C02, C;^i, C20), CTq = (/ill, ^20' ^4o)' 

where 60 and ctq satisfy 

(^20)^ + ^40 ^ 0' ^00 ~ ~^02' ^02 7^ 0; 
^20 ~ ^11^*1 ~ 2Co2[(/iii)^ — /i2o]- 

Then, for some [e, 5, a) near (0, ^O; o"o) ^/^e system fl7T|) /ias 2 /zmzt cycles near the origin. 

Proof. Let M(/;,, 5) be the Melnikov function of system flTTl) near the origin. Then, by 
Theorem 2, we have 

M(/i,(5) = /i3/^^6,((5)/i^/l 

i>o 

Executing the Mathematica program yields 



_ 2r(5/4)V27r 

^° - r(7/4)(5- 2^2)1/4^"°° + 

Letting cqo = — C02 (implying 60 = 0), we obtain the following expressions for 6^'s: 

y'z^^ -^/tt 

feiko=o = r(9/4)(-2/.io -"2/^40)3/^^''° ~ + ~ 

I _ r(l/4)0r \ X 1 ('721 

&2ko=o - 4r(ii/4)(-2/.io- 2/^40)^/^^ ''^"'° " ""''^ ^ 

63160=0 = r(7/4)v/7r [53i(c2o - /inCn) + ^32020], 



where 





3hi 
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12hl,h2o + 4:hlf^ - 8h. 



'32 



3/1? 



11 



15/ili/i2o + 8/i?i(/i^o - 2/i4o) + 12/^^0 + 24/i2o/i. 



■40 • 



Let 



bj = bj\bo=o, j = 1,2,3. 



Noticing (172]) . we can obtain C20 = huCu — 2co2(/iii — /i2o) by solving 61 = 0. Substituting 
C20 = hiiCii — 2co2(/?.ii — /i2o) into 62 and 63 results in 



For (/120, ^4o) near (/i205^4o)' have —2/^20 ~ 2/140 > 0. Then 62 = if and only if 
C02 = 0, which yields 63 = 0. Hence, it is clear that the conclusion follows from Theorem 



The proof is complete. 



Therefore, system (l64l) has 2 limit cycles near the nilpotent center (0, 1). For the sym- 
metry, system ( IMIl also has 2 limit cycles near the nilpotent center (0, —1). Namely, there 
are altogether 4 limit cycles of system f l6^ bifurcating from both symmetric nilpotent 
centers (0, ±1). 



Appendix 

Basing on the formulae given in section 2 we have written the Mathematica code for 
computing {bj}, presented below. It contains several subroutines (as shown in the code) 
for computing ej, qj{x), H*{x), aj{x), qj{x), ipi^x), qj, Vij, (3ij and bj. 



########## compute the e_j=e[j] coefficients ########## 

H=h [2 , 0] x-2+h [0 , 2] y-2+Sum [Sum [h [i- j , j ] x-{i- j }y- j , { j , , i>] , { i , 3 , n+1}] ; 

div=Sum [Sum [c [i-j , j] x"{i-j}y~ j ,{j ,0, i}] ,{i,0,n-l>] ; 

phi=Sum[e [j]x~j ,{j ,2,m}] ; 

H1=D [H , y] ; H1=H1/ . ■[y->phi} ; 

Do [H2[i]=Coefficieiit [Hl,x,i] ,{i,2,m}] ; 

Do [t [i] =Solve [H2 [i] ==0 , e [i] ] , {i , 2 ,m}] ; 

?t 

########## compute q_j(x)=q[j] ########## 

Do [q [ j ] = ( 1/ j ! ) D [div , {y , j - 1}] , { j , 1 , n}] ; 

Do [q [ j ] =q [ j ] / . {y->plii} , { j , 1 , n}] 

########## compute H*_j (x)=Hs [j] ########## 

Do [Hs [j]=(l/(j+l) !)D[div,{y,j+l}] ,{j,l,n}] ; 

Do [Hs [ j ] =Hs [ j ] / . {y->phi} , { j , 1 , n}] 

########## compute a_j(x)=a[j] ########## 

Vl=Sum[a[j]w-j ,{j ,l,m}] 
Ht=Sum[Hs [j]V-{j+l},{j ,1,11}] 
F=Ht-w~2; 

Do [W[j]=Coeff icient [F,w, j] ,{j ,l,m}] 

Do [tl[j]=Solve [W[j+l]==0,a[j]] ,{j,l,m}] 

?tl 

########## compute bar_q_j (x)=qb [j] ########## 
Q=Sum[q[j]Vl-j,{j,l,ii}] 



r(i/4)v/^ 



r(ll/4)(-2/.2^- 2/^40)1/^ 
-16r(7/4)v/^ {hi, - 2/l2o)(/i^o + ^40)C02. 



3. 



Do [qb [j ] =2Coef f icient [Q , w , 2 j + 1] , { j , ,m}] 
### compute psi(x)=psi ### 
HO=H/.{y->phi} 
psi=[H0]-{l/(2p)}; 

########## compute tilde_q_j=qt[j] 
Psi=Sum [PS [ j ] u" j , { j , 1 , m}] 
HOO=HO/.{x->Psi} 
G=H00-{2p}-u 

Do[g[j]=Coefficient[G,u,j] ,{j ,0,m>] 
Do [t2 [ j ] =Solve [g [ j ] ==0 , PS [j ] ] , { j , 1 ,m}] 
?t2 

Do[qt[j]=(qb[j]/.{x->Psi})D[Psi,u] ,{j ,0,m}] 

########## compute r_ij coefficients ########## 

Do[r[i,j]=2Coeff icient [qt [j] ,2i] ,{i,0,m},{j ,0,m}] 

########## compute beta_ij=bt [i, j] coefficients ########## 

Do [bt [i, j]=Integrate [u~{(l-p+2i)/p>(l-u"2) "j Sqrt [l-u~2] ,{u,0,l}] ,{i,0,m},{j ,0,m}] 

########## compute b_j coefficients ########## 

Do [b [ j ] =Sum [r [ j -p*i , i] bt [ j -p*i , i] , {i , , IntegerPart [ j /p] }] , { j , , m}] 



Executing the above program yields expressions in the original coefficients hij and Cij. 
Below we list the final coefficients bj (the intermediate expressions such etc. are 

omitted) for p = 2, a; = 1 and n = 3. We found 



2tt 



(4/i40-2/iii)V4 



Coo- 



If bo — 0, i.e., Coo — 0, then 
16 



bi 



3 (4/i4o - 2/iii)7/4 
-/l2l(4/l4o - 2/i^Jcoi 



(4/i4o - 2/iii)c20 - (/il2^21 - 3/i21^31 + 4/ii2/i4o)ClO 



62 = 



TT 



X 



|2(/i|i - 2/i4o) - 12/i|i/l22 + 30/ii2/i2l/i3l(/iii + 4/i40) 



(4/140 - 2/li,)lV4 - rV-21 
+/l^l(-35/lii + 8/l22/i40) + 4/l40(-5/lii + 8/i22/i40) + 3/i?2(^21 " 24/l2^/i4o - 16/i|o) 



-4/103/^21 (/ill - 14/12^/140 + 2Ahlo) 

— 32(/l2i — 2/l4o)^/l4oCo2 



C20 



-8(/iii - 2/i40)^ /il2/i2l(/i21 - 12/i40) + /i3l(3/iii + 4/i4o) 



Cll 



+2(/i2^ - 2/i4o) [4/io3(/iii - 2/i4o)(/iti - 24/ilo) + 6/112/131 (/^si - 24/iii/i4o - 16/ilo) 

+5/l?2^2l(/i21 - 8/l^i/l40 + 48/l^o) + /i2l(-4/l^i/l22 + 12/l4o(5/lli - 8/l22/i40) 



+/i2i(15/l|i + 56/l22/i4o)) 



Coi 



-8/113/121 (/ill - 12/i4o) (/ill - 2/i4o)^ - 15/1^2/121/131 (/i2i - 40/iii/i4o - 80/1^0) 



-5/if2(/i2i - 12/i|i/i4o + UAhl^h'lQ + 64/i|o) + 3/i3i(4/io3(/i2i - 26/ifi/i4o + 32/i2^/i: 



40 



+32hlo) + 5/i2i(4/iti/i22 + 4/i4o(3/iii - 8/122/^40) + ^ii(7/iii + 8/122/^40))) 

+/ll2(20/lo3/i2l(/lii - 2/l40)(^21 - 8/lii/l40 + 48/l|o) + 3(4/i^i/l22 + 16/i2o(-5/i 
+8/i22/l40) + 8/lii/i40(-35/lii + 16/i22/i40) " /i2l(35/iii + 104/l22/i40)) Cio}, 

16 



15(4/i40-2/i2 )27/4 



x(-..) 
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